Anisotropic behaviors of massless Dirac fermions in graphene under periodic potential 
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Charge carriers of graphene show neutrino-like 
linear energy dispersions as well as chiral behav- 
ior near the Dirac pointi*2i2i£i&i&. Here we re- 
port highly unusual and unexpected behaviors of 
these carriers in applied external periodic poten- 
tials, i.e., in graphene superlattices. The group 
velocity renormalizes highly anisotropically even 
to a degree that it is not changed at all for 
states with wavevector in one direction but is re- 
duced to zero in another, implying the possibil- 
ity that one can make nanoscale electronic cir- 
cuits out of graphene not by cutting it 7,8 ' 9 ' 10 but 
by drawing on it in a non-destructive way. Also, 
the type of charge carrier species (e.g. electron, 
hole or open orbit) and their density of states 
vary drastically with the Fermi energy, enabling 
one to tune the Fermi surface-dominant proper- 
ties significantly with gate voltage. These re- 
sults address the fundamental question of how 
chiral massless Dirac fermions propagate in peri- 
odic potentials and point to a new possible path 
for nanoscale electronics. 

Since the pioneering work by Esaki and Tsui 1 -, super- 
lattices have been studied extensively and have had a 
huge impact on semiconductor physic o 12 ' 13 . Superlattices 
demonstrate interesting phenomena such as negative dif- 
ferential conductivity, Bloch oscillations, gap openings 
at the mini Brillouin zone boundary formed by the addi- 
tional periodic potential, et o 12 ' 13 . Conventional semicon- 
ducting superlattices are mainly produced by molecular- 
beam epitaxy and metallo-organic chemical vapour-phase 
deposition while metallic superlattices are made by sput- 
tering procedure o 12 ' 13 . We expect that, by modulat- 
ing the potential seen by the electrons, graphene su- 
perlattices may be fabricated by adsorbing adatoms on 
graphene surface through similar techniques, by position- 
ing and aligning impurities with scanning tunneling mi- 
croscop y 14 ' 15 i 16 , or by applying a local top-gate voltage 
to graphen e - 17 ' 18 ! 19 . Epitaxial growth of graphene^ on 
top of pre-patterned substrate is also a possible route to 
graphene superlattice. Recently, periodic pattern in the 
scanning tunneling microscope image has been demon- 
strated on a graphene monolayer on top of a metallic 
substrat o 20 i 21 ' 22 as well. 

The low energy charge carriers in pristine graphene are 
described by a massless Dirac equation and have a linear 
energy dispersion which is isotropic near the Dirac points 
K and K' in the Brillouin zon(£&&2&21&. (Fig. la). It is 
shown experimentally that the carriers have a group ve- 



locity of vo 10 6 m/s which plays the role of an effective 
speed of light in (2+1) dimensional quantum electrody- 
namics 3,4 . Within the effective-Hamiltonian approxima- 
tion, 



H = Hvq 
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where k x and k y are the x and y components of the 
wavevector k of the Bloch state defined with respect 
to the Dirac point, the wavefunction of the quasiparti- 
cles in graphene has two components corresponding to 
the amplitude on the two different trigonal sublattices of 
graphene and can be represented by a two component 
spinor 23 ' 24 ' 25 . This spinor structure of the wavefunction 
is called a pseudospin (because it is not related to a real 
spin) or chiralit y 23 ' 24 ' 25 ' 26 ' 27 , which is of central impor- 
tance to the novel physical properties of graphene super- 
lattices discussed below. 

Let us now consider the situation that an additional 
periodic potential is applied to graphene. If the spatial 
period of the superlattice potential is much larger than 
the nearest neighbor carbon-carbon distance in graphene 
(~ 1.42 A), the scattering of a state close to one Dirac 
point to another one does not occu r 23 i 24 i 26 . Therefore, 
even though there are two nonequivalent Dirac cones for 
the energy dispersion surface of graphene, focusing on a 
single cone is sufficient. This condition also implies that, 
in the graphene superlattices discussed here, there is no 
gap opening at the Dirac poin t 23 ' 24 ' 26 . 

To investigate the physics of charge carriers in 
graphene superlattices, we have calculated the en- 
ergy dispersions, the group velocities, and energy gap 
openings at the minizone boundaries (MB) within the 
effective-Hamiltonian formalism 1 -. Effects of the external 
periodic potential are incorporated into our calculation 
through the scattering matrix elements between pseu- 
dospin states, or chiral eigenstates, of the electrons in 
graphene: 



( S ,k|[/(r)| S ',k')=^i(l 
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Here, G's are reciprocal lattice vectors of the superlattice 
and #k,k-G is the angle between k and k — G. U(r) and 
U(G) are external potential in real and wavevector space, 
respectively, s and s' are either +1 or —1 depending 
on whether the energy of the state is above or below 
the energy at the Dirac point, respectively. We have 
also carried out a tight-binding formulation and obtained 
identical results as those discussed below. 
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FIG. 1: Graphene superlattices and anisotropic Dirac 
cones, a, Schematic diagram of graphene. Inset: the Bril- 
louin zone of graphene and Dirac cones centered at Dirac 
points among which two (K and K') are nonequivalent (left) 
and the linear and isotropic energy dispersion near one of the 
Dirac points of charge carriers in graphene (right), b, One- 
dimensional (ID) graphene superlattice formed by Kronig- 
Penney type of potential periodic along x direction with spa- 
tial period L and barrier width w. The potential is Uxd in 
the grey regions and zero outside. Inset: energy dispersion 
of charge carriers in ID graphene superlattice. The energy 
dispersion along any line in two-dimensional (2D) wavevector 
space from the Dirac point is linear but with different group 
velocity. For a particle moving parallel to the periodic direc- 
tion, the group velocity (v\\ ) is not renormalized at all whereas 
that for a particle moving perpendicular to the periodic di- 
rection (v±) it is reduced most, c, 2D graphene superlattice 
with muffin-tin type of potential periodic along both x and y 
directions with spatial periods L x and L y , respectively. The 
potential is U20 inside the grey disks with diameter d and 
zero outside. Inset: energy dispersions of charge carriers in 
2D graphene superlattice. 



First, for a one-dimensional (ID) graphene superlattice 
(Fig. lb), we find that the group velocity for states with 
wavevector k (k is the wavevector of the Bloch state de- 
fined with respect to the Dirac point), is anisotropically 
renormalized, i.e., it is a strong function of the direction 
of k. For pristine graphene, the group velocity of states 
near the Dirac point is parallel to k and of constant mag- 
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FIG. 2: Anisotropic velocity renormalization in 
graphene superlattices. a, The component of the group 
velocity parallel to the k vector [uj = v(k) ■ k with k mea- 
sured from the Dirac point] of charge carriers in a ID graphene 
superlattice in units of the Fermi velocity in graphene (vo) 
versus the angle (#k) of the k- vector from the periodic poten- 
tial direction x, (solid lines) and that in a superlattice made 
from a fictitious system of nonchiral fermions with proper- 
ties otherwise identical to those in graphene (dashed lines). 
Red, green and blue lines correspond to Uro being 0.2 eV, 
0.3 eV and 0.5 eV, respectively, b, Similar quantities as in a 
for a rectangular 2D graphene superlattice. Red, green and 
blue lines correspond to L^d being 0.3 eV, 0.5 eV and 0.7 eV, 
respectively, c, The group velocity of charge carriers in a 
ID graphene superlattice (solid line) with k perpendicular to 
the periodic direction, v_i_, in units of vo versus [/id (solid 
line) and that in a superlattice made from a fictitious sys- 
tem of nonchiral fermions with properties otherwise identical 
to those in graphene (dashed line), d, v± versus the poten- 
tial spatial period (L) of charge carriers in a ID graphene 
superlattice. Red, green and blue lines correspond to a fixed 
potential barrier height but with width (w) being 5 nm, 10 nm 
and 25 nm, respectively. 



nitudc (vq). For example, in a ID superlattice of Kronig- 
Penny type of periodic potential with potential barrier 
height (E/id) of 0.5 eV and spatial period (L) and barrier 
width (w) of 10 nm and 5 nm, respectively, the group 
velocity of the charge carriers when k is along certain 
direction is renormalized to be slower than 40 % of its 
original value vq but is the same as vq along some other 
direction. [Fig. 2a: the plotted quantity Vt = v(k) • k 
is the component of the group velocity parallel to the 
wavevector k in units of Vq. We note that this quantity 
which depends only on the direction of k (Supplementary 
Discussion 2) is exactly equal to the absolute value of the 
group velocity v g when k is at 0, 90, 180 or 270 degrees 



from the periodic direction of the applied potential and, 
when the applied potential is weak, is only slightly differ- 
ent from v g at other angles (Supplementary Discussion 
3).] Thus, the group velocity of charge carriers can be 
tailored highly anisotropically in graphcne supcrlattices. 
More interestingly, the group velocity when k is along 
the direction perpendicular to the periodic direction of 
the potential (v±) is reduced the most, whereas when k 
is in the parallel direction, it is not renormalized at all 
(Fig. lb) . This result is counter- intuitive since the veloc- 
ity is strongly reduced when the charge carrier is moving 
parallel to the hurdles, but is not modulated when it is 
crossing them. 

To understand the physics behind this phenomenon, 
we have performed the same calculation for a fictitious 
system with carriers that have no chirality but other- 
wise identical to those in graphene including the linear 
energy dispersion. The group velocity in this system is 
reduced isotropically and the renormalized group veloc- 
ity is close to v±, i.e., the maximally renormalized one, 
in ID graphene superlattices (Fig. 2a). Thus, it is clear 
that the absence of velocity renormalization in the direc- 
tion parallel to the periodicity of the external potential 
originates from the chiral nature of the electronic states 
of graphene. This behavior can be demonstrated more 
clearly by second order perturbation theory in the case 
of the ID periodic potential with weak amplitudes (Sup- 
plementary Discussion 2). We note that the chirality dis- 
cussed here is also important in tunneling phenomenon 
in graphene through a single barrier— or a finite number 
of barriers 2 ^. 

In the case of two-dimensional (2D) graphene super- 
lattices, the group velocity is renormalized for k states 
along every direction, but anisotropically (Fig. lc). As 
the amplitude of the potential increases the overall group 
velocity is reduced and the ratio of the maximum group 
velocity to the minimum one is enhanced (Fig. 2b). Here, 
again, the anisotropy disappears if the chiral nature of 
the states in graphene is arbitrarily removed. As was 
demonstrated for the ID superlattice, the sinusoidal de- 
pendence on the angle of propagation as well as the over- 
all shift in the case of 2D graphene superlattice of the 
component of the renormalized group velocity parallel to 
k in the weak potential limit is well explained and repro- 
duced by second order perturbation theory (Supplemen- 
tary Discussion 2). 

Remarkably, the anisotropy in energy dispersions of 
the ID superlattices can be tuned by changing the ap- 
plied potential in such a way that v± is reduced com- 
pletely to zero (Fig. 2c). Hence, we can achieve ex- 
tremely low mobility in one direction and normal con- 
duction in another one simultaneously. This enables one 
to control the flow of electrons dramatically. It also pro- 
vides a novel non-destructive pathway to make graphene 
nanoribbons^&^Jii which have been actively pursued by 
way of cutting graphene sheets*^. The chiral nature 
of the states in graphene also plays a decisive role here. 
In the model without chirality as discussed before the 
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FIG. 3: Energy gap at the superlattice Brillouin zone 
or minizone boundary of a ID graphene superlattice. 

a, Energy of charge carriers in ID graphene superlattice ver- 
sus the component of the wavevector k parallel to the periodic 
potential direction (k m ) at a fixed k y . Dashed vertical lines 
indicate minizone boundaries (k x = ±7r/L). AE is the en- 
ergy gap at the minizone boundary for a given ky. Red and 
blue lines correspond to k y being zero and 0.012 A" 1 , respec- 
tively, b, AE versus k y for charge carriers in ID graphene 
superlattice (solid lines) and that in a superlattice made from 
a fictitious system with states without chirality but otherwise 
identical to graphene (dashed lines). Red, green and blue 
fines correspond to Uio being 0.1 eV, 0.3 eV and 0.5 eV, 
respectively. 

(isotropic) group velocity of charge carriers is reduced 
monotonically and never reaches zero within a conceiv- 
able range of the potential amplitude (Fig. 2c). We can 
also achieve vanishing group velocity in one direction by 
changing the length parameters of the superlattice (Fig. 
2d). 

Graphene superlattices show peculiar behavior of gap 
openings at the MB formed by the external periodic po- 
tential (Fig. 3). In conventional layer-structured ID su- 
perlattices, gap opening at the MBs is considered to be 
nearly constant, independent of k. ID graphene super- 
lattices, however, are different in that the gap (AE) van- 
ishes when k is along the direction of the periodic po- 
tential, i.e., at the centre of the MB (Fig. 3a and 4a). 
Moreover, the size of the gap depends strongly on where 
it is on the MB (Fig. 3b). These strong anisotropics in the 
gap opening do not happen in superlattices made from a 
system having linear energy dispersions but no chirality 
(Fig. 3b). Hence, again, the chiral nature of charge car- 
riers in graphene is key in generating these anisotropics 
in the gap opening as it does in the velocity renormaliza- 
tion. In particular, the gap closure at the centre of the 
MB is directly related to the absence of back-scattering 
of charge carriers from a scattering potential when the 
size of the scatterer is several times larger than the inter- 
carbon distanc o 23 ' 24 i 26 . In ID graphene superlattices, the 
important length-scale is L, which is much larger than the 
inter-carbon distance, and hence the gap does not open 
at the centre of the MB. 

The largest gap at the MB in a graphene superlattice 
is proportional to the amplitude of the applied potential 
if the potential is weak (i.e., small compared to the band 
width) and the size of which thus can be made to be a few 
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FIG. 4: Energy dispersions and densities of states of 
charge carriers in graphene superlattices. a, Energy 
of charge carriers in ID graphene superlattice with Uio = 
0.3 eV, L = 10 nm and w = 5 nm in the first (red and black) 
and the second (blue and pink) band above the vertex of the 
Dirac cone versus 2D wavevector k away from the Dirac point. 
Minizone boundaries are at k x = ±0.031 A -1 . Arrows indi- 
cate points on the minizone boundary where the gap closes, 
b, Density of states (DOS) of charge carriers in electron or- 
bits (red), open orbits (green) and hole orbits (blue) in the ID 
graphene superlattice characterized in a versus the Fermi en- 
ergy (E-p). The origin of the energy scale is set at the energy 
of the Dirac point. The DOS of each species is the height of 
the corresponding colored region. Dashed black line shows the 
DOS of pristine graphene for comparison, c, Similar quanti- 
ties as in a for a 2D graphene superlattice with U20 = 0.3 eV, 
L x = L y — 10 nm and d = 5 nm. d, Similar quantities as in 
b for the 2D graphene superlattice specified in c. 

tenths of an electron volt with appropriate perturbation 
and much larger than room temperature (Fig. 3b, Sup- 
plementary Discussion 5) . We have also investigated the 
gap opening in ID graphene superlattices with different 
values of length parameters (L and w). We find that, by 
changing these parameters, the anisotropy in the gap at 
the MB can also be controlled (Supplementary Discus- 
sion 5). 

Due to the velocity renormalization near the Dirac 
point and the strong anisotropy in energy dispersion close 
to the MB, the type and the density of states (DOS) of 



charge carriers vary drastically from those in graphene as 
one varies the Fermi energy (Fig. 4). As the Fermi level 
changes, the topology of the Fermi surface also exhibits 
a dramatic variation (Supplementary Discussion 6). For 
example, as the Fermi level increases from the energy at 
the Dirac point, the charge carriers of a ID graphene su- 
perlattice fill electron orbits and show a linear increase 
in the DOS with slope larger than that of graphene; 
but above certain value, the DOS of electron orbits van- 
ishes and charge carriers suddenly fill open orbits and 
hole pockets. When the Fermi level increases further, 
charge carriers are in purely open orbits and then the 
DOS of electron orbits starts to reappear and increases 
again (Fig. 4b). We expect that the Fermi level in a 
graphene superlattice can be tuned as in graphene by 
applying a gate voltag o 3 i 4 i 18 ' 19 . Hence, by exploiting the 
various characteristics of charge carriers and the Fermi 
surface topology, one can manipulate a variety of phys- 
ical properties dominated by the Fermi surface, such as 
conductance or magnetoresistance, significantly. 

The anisotropic gap opening at the MB and the dra- 
matic variation of the characters of charge carriers with 
the Fermi energy are also common in 2D graphene super- 
lattices. The gap at the centres of the zone boundaries 
closes as in ID graphene superlattice (Fig. 4c). However, 
the gap at the corners of the 2D MB also disappears. 
This behavior, which occurs in rectangular 2D graphene 
superlattices in general, has again its origin in the chi- 
ral nature of charge carriers in graphene (Supplemen- 
tary Discussion 4). In a square 2D graphene superlattice, 
charge carriers are electrons, holes or a mixture of the two 
depending on the Fermi level (Fig. 4d). For general rect- 
angular 2D graphene superlattices, charge carriers can 
also be in open orbits. 

Here we have presented several novel physical proper- 
ties of graphene superlattices with Kronig-Penney type 
ID and muffin-tin type 2D potentials. Through addi- 
tional calculations, we have confirmed that all the salient 
features of our findings are the same in sinusoidal or 
Gaussian types of graphene superlattices in general as 
well. The novel properties discovered in the present study 
thus should be of interest to the fundamental study and 
practical applications of graphene systems in general. 

Finally, since the massless Dirac fermions in graphene 
superlattices have some features in common with high- 
energy relativistic particles propagating in anisotropic 
space like the anisotropy in the group velocity^, inter- 
esting physics of the latter may also be investigated by 
table-top experiments based on our theoretical findings. 
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Supplementary Discussion 1 : Effective-Hamiltonian 
formalism 



energy zero, respectively. Eigenstates of this Hamiltonian 
is given by 



There are two carbon atoms per unit cell in graphene, 
forming two different sublattices, and hence the eigen- 
state of charge carriers in graphene can be represented 
by a two component basis vector. The Brillouin zone of 
graphene is a hexagon which has two inequivalent ver- 
tices, so-called the Dirac points, K and K', that cannot 
be connected by reciprocal lattice vectors. In this work, 
we are considering eigenstates near K only as discussed 
in the paper. The effective Hamiltonian for low-energy 
quasiparticles of graphene in this basis is given by 



#o(k) = hv 







o 



(i) 



where vq is the Fermi velocity and k the small wavevector 
of the quasiparticle from the K point in the Brillouin zone 
of graphene. The energy spectrum of this Hamiltonian is 
E = shvak where s is +1 or —1 for an eigenstate above or 
below the Dirac point energy which is defined to be the 



(r| S ,k) 



1 



ol (K+k)-r 



1 

ise lt 



(2) 



where #k is the angle of vector k with respect to the k x 
direction. Now, when an additional periodic potential 
U (r) is applied to graphene, the scattering amplitude 
between states are given by 

(a, k \U(t)\ s', k') = \ (1 + ss'e-^-*) U(G) <5 k ,, k _ G , 

G 

(3) 

where G and U(G) are the reciprocal lattice vector and 
the corresponding Fourier component of the external pe- 
riodic potential, respectively, and 0k,k-G the angle from 
k G to k. The energy dispersions and eigenstates of 
the quasiparticles in a graphene superlattice are obtained 
non-perturbatively within the single-particle picture by 
solving the following set of linear equations: 



J 



(E - e s>k ) c(s, k) = \ I 1 + ss'e-^-v) U(G) c{s', k G) 



(4) 



s',G 



r 



where E is the superlattice energy eigenvalue and e Si k — 
shvok the energy of the quasiparticles before applying the 
periodic potential. c(s, k) and c(s', k— G) are the ampli- 
tudes of mixing among different unperturbed quasiparti- 
cle states. 



Supplementary Discussion 2 : Velocity 
renormalization near the Dirac point from second 
order perturbation theory 

When the external potential is weak, perturbative cal- 
culations can give results in good agreement with the full 
calculation and also in physically more intuitive forms. 
For pristine graphene, the group velocity of states near 
the Dirac point is parallel to k and of constant magnitude 
(i>o). For a graphene superlattice, the renormalization in 
the component of the group velocity of quasiparticles par- 
allel to the wavevector k [vt = v(k) • k] around the Dirac 
point can be obtained within second order perturbation 
approximation as 



f'O 



G#0 



«g|G|» 



'k.G 



(5) 



where #k,G is the angle from G to k. From Eq. (T5J), it 
is clear that for weak potentials, the velocity renormal- 



ization grows as square of the amplitude of the exter- 
nal potential and Vf, depends only on the direction of k. 
Throughout this document, by weak potential we mean 
that the condition, 



\U(G)\ 
v Q \G\ 



«1 , 



(6) 



is satisfied for all the U(G) components. In the absence 
of chirality of the states, the scattering matrix element 
[Eq. ©] should be changed into 



(s,k|£/(r)| S ',k')=^(G) <W- G 



(7) 



G 



Using Eq. ([7]), the similar quantity as in Eq. (JSJ) for 
nonchiral massless Dirac fermions is now given by 



non— chiral 



£ 

G^O 



^|GP 



(8) 



which is isotropic independent of the direction of k. Com- 
paring Eq. ([5]) and Eq. ([8]), the renormalization of group 
velocity in a one-dimensional (ID) superlattice for a ficti- 
tious graphene with states without chirality corresponds 
to the maximum renormalization in the corresponding 
ID graphene superlattice. This trend agrees with the re- 
sults from the full calculation when the potential is weak 
(Fig. 2 of paper). 
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FIG. 5: (Supplementary Figure 1) Dependence of the 
velocity renormalization on the amplitude of periodic 
potential in a ID graphene superlattice. Square root of 
the difference between the group velocity for state with k 
along the direction perpendicular to the periodic direction of 
the potential (v±) and the unrenormalized one (vq) divided 
by vo versus the potential amplitude Uiu- Solid red line and 
dashed blue line are results of the full calculation and second 
order perturbation theory, respectively. 



If Eq. ([5]) is applied to the Kronig-Penney type of ID 
graphene superlattice periodic along the x direction as 
discussed in the paper (Fig. lb of paper), 



Vo 



El . 2 f^W 



n>0 



7k, x 



(9) 



where L is the spatial period of the potential, and Urn 
and w are the height and the width of the rectangular 
potential barrier, respectively. Equation (|9|) is in good 
agreement with the full calculation when the potential is 
weak (Supplementary Fig. 1) and also shows sinusoidal 
behavior with respect to the polar angle of k with respect 
to the periodic direction of the potential as well as the 
absence of renormalization for a particle moving across 
the potential barriers. 



For the case of muffin-tin type of two-dimensional (2D) 
' ice (Fig. lc of paper) periodic along 



graphene superlattice (Fig. 

' i: * : : * •• -~ods L x and L y , respec 



-_j — i/VHOOUl^ ig,. 1^ ^-L j 

both x and y directions, with perio 
tively, 



Vo 



vo 



2ir 2 U 2 D d 2 

„2r2r2 



sr 1 T *( Gd 

2^ G 4 J l \ 9 



sin 6> k ,G 



(10) 

where U2D is the height of the circular potential barrier 

with diameter d and G = f ff-w, j^-n^j the reciprocal 

lattice vector with integers m and n and J\ (x) is the first 
order Bessel function of the first kind. When L x « L y 
(as in Fig. 2b of paper), the dominant contribution in 
the sum comes from the terms with |m| = l,n = and 
m = 0, \n\ = 1. Counting only these four terms, 



"0 




k , x 



(ii) 



where the relation sin 2 #k,y = 1 — sin 2 0k, & has been used 
in the second line. Now the group velocity v^. is renormal- 
ized in every direction. Equation (jTTJ) reproduces the si- 
nusoidal variation of the velocity renormalization as well 
as the constant shift as obtained in the full calculation 
quite well when the potential is weak (Fig. 2b of paper). 



value of the group velocity v g when k is at 0, 90, 180 
or 270 degrees from the periodic direction of the applied 
potential and, when the applied potential is weak, is only 
slightly different from v g at other angles (Supplementary 
Fig. 2). 



Supplementary Discussion 3 : The magnitude and 
the component parallel to the wavevector k of the 
group velocity 

The component of the group velocity parallel to 
the wavevector k above is exactly equal to the absolute 



Supplementary Discussion 4 : Band gap at the 
minizone boundary from degenerate perturbation 
theory 

When the wavevector k is on the minizone boundary 
(MB) of the ID graphene superlattice, k x = ±tt/L, two 
states |s,k) and |s,k— (2tt/L,0)) are degenerate before 
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applying the periodic potential. The largest contribution 
to the energy eigenvalues at the MB comes from these two 
degenerate states. Scattering amplitude between these 
two states is given by Eq. ([3]). By choosing the origin 
at the centre of a potential barrier, we can make all the 
Fourier components of the periodic potential real without 
losing generality. Now the Hamiltonian for the two-state 
system is 



FIG. 6: (Supplementary Figure 2) The magnitude and 
the component parallel to the wavevector k of the 
renormalized group velocity in a 2D graphene super- 
lattice. The component of the group velocity parallel to the 
k vector [uj, = v(k) • k with k measured from the Dirac point] 
of charge carriers in a 2D graphene superlattice (solid lines) 
and the absolute value of the group velocity (dashed lines) 
in units of the Fermi velocity in graphene (vq) versus the an- 
gle (#k) of the k-vector from the periodic potential direction 
x. Red, green and blue lines correspond to the potential am- 
plitude C/2D being 0.3 eV, 0.5 eV and 0.7 eV, respectively. 
Plotted quantities are obtained from the full calculation by 
solving Eq. (J4j . 



J 



H(k) 



i ( 1 + e '» k , M 2,/i, )) U(2w/L) 



l + e-'fW^/L.o)) U(2ir/L) 

£ s,k-(27r/L,0) 



(12) 



where £ Si k = £ S! k-(27r/L,o) = shv^k is the energy of the 
charge carrier before the external potential is applied and 
U (2ir/L) the Fourier component of the periodic potential 
whose wavevector connects the two MBs at k x — ±ir/L. 
The energy separate of the eigenvalues of Eq. (fT2"|) (i.e., 
the energy gap) is given by 



AE = 2 \U(2ir/L) sin 9 K£ \ 



(13) 



where #k, s is the polar angle between k and x. Equa- 
tion (p~3]) clearly shows that the gap opening depends 
on k and, in particular, is zero at the centre of the MB 
(Fig. 2 of paper), and that, as discussed in the paper, the 
maximum gap opening is proportional to the amplitude 
of the external potential in the weak potential regime. 

If the chirality of the states in graphene is manually re- 
moved by using Eq. ([7]) for the scattering matrix element, 
the energy gap becomes 



A_E Ilnn _ 



chiral 



2 \U(2it/L)\ 



(14) 



in which case the gap neither closes at the centre of the 
MB nor depends sensitively on the position along the MB 
in the weak potential limit (Fig. 2 of paper). 

For a 2D rectangular graphene superlattice, in 
which the primitive translational lattice vectors are 
orthogonal, for the same reason as in the ID case, 



the gap closes at the centres of MBs, i.e., when k is 
at (±tt/L x ,0) or (0,±7r/Lj,). The more interesting 
case is the corner of the minizone, where the four 
degenerate states \s,(n/L x ,Tr/L y )), \s, (— ir/L x , ir/Ly)), 
\s,(tt/L Xi ~tt/L v )) and \s,(-Tr/L x ,-Tr/L y )) mix 
strongly among themselves by the applied periodic 
potential. If we set up a similar matrix for this case like 
Eq. (|12[) for the ID graphene superlattice, the energy 
eigenvalues of the matrix are 



E : 



± 



[U(2n/L x ,0) sin 9 Kx f - [£7(0, 2 7 r/L J/ )cos 6 Kx ] 

(15) 



where k = (n / L Xl ir / L y ) is at the minizone corner and 
each eigenvalue is doubly degenerate. The energy spec- 
trum given by Eq. (|15p clearly shows that there is no 
gap at the minizone corner between the first and the 
second band (Fig. 4c of paper). This gap closure at 
the minizone corner is not obvious because a transition, 
which is not of a backscattering process, from one of the 
four k points at the zone corners to another can oc- 
cur in the 2D rectangular graphene superlattice. For 
example, the state |s, (ir/L x , 7r/L y )} can mix with an- 
other state |s, (— tt/L x , n/L y )) by the reciprocal lattice 
vector G = (2ir/L x ,0) but the two k vectors are not 
anti-parallel. To understand the origin of the gap closure 
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at the minizone corner, we repeated a similar calculation 
for a 2D rectangular superlattice formed of a non-chiral 



J 



system with linear energy dispersions, 
energy eigenvalues are 



In this case, the 



non— chiral 



e s , k + U(2tt/L x , 2n/L y ) ± [U(2n/L x , 0) - [7(0, 2ir/L y )] 
e s , k - U(2tt/L x , 2ir/L y ) ± [U(2w/L x , 0) + £7(0, 2ir/L v )] 

I 



(16) 



According to Eq. (fTfj)) . there is a finite energy gap be- 
tween the first and the second band in general, other than 
in accidentally symmetric cases. Therefore, the gap clo- 
sure at the minizone corner of a 2D rectangular graphene 
superlattice is a direct consequence of the chiral nature 
of the states in graphene. 



Supplementary Discussion 5 : Dependence of the 
band gap at the minizone boundary on length 
parameters and broken particle-hole symmetry 

For a Kronig-Penney type rectangular potential bar- 
rier ID graphene superlattice, the energy gap at the MB 
can be expressed with Eq. (fl"3)) in terms of superlattice 
parameters as 



AE= - 

7T 



■I 7TW 

UiB sin ( — ] sin y k ,i 



(17) 



Thus, according to this simple degenerate perturbation 
theory result if the spatial period (L) becomes long for 
a constant barrier width (w) the gap should decrease, 
whereas if L is fixed and w is increased from zero, the 
gap should reach a maximum at w = L/2 and after that 
should decrease symmetrically. The former trend is ob- 
served in the full calculation (Supplementary Fig. 3a); 
however, the latter seems not to hold in the full cal- 
culation (Supplementary Fig. 3b). Moreover, the gap 
openings at the MB above and below the Dirac cone are 
different (Supplementary Fig. 3), which shows the limi- 
tation of the simple degenerate perturbation theory be- 
cause the potential is not so weak and possesses higher 
Fourier components. One thing to note is that the energy 
dispersion for states above the energy of the Dirac point, 
including the gap opening, for a ID Kronig-Penney type 
superlattice with width w = wq is identical with that 
of states below the energy of the Dirac point for width 
w = L — wq (Supplementary Fig. 3b, compare red and 
blue lines). This symmetry can be understood follow- 
ing a simple argument. If we start from a ID Kronig- 
Penney superlattice with w = wq and then change w to 
be w = L — wq and at the same time invert the whole 
potential, the resulting periodic potential is identical to 



the original one other than a constant shift which may 
be ignored. Inverting the potential is effectively the same 
as exchanging particles with holes. 
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FIG. 7: (Supplementary Figure 3) Dependence of the 
energy gap at the minizone boundary on the length 
parameters of ID graphene superlattice. a, The energy 
gap AE between the first and the second band at the minizone 
boundary versus k y for charge carriers above (solid lines) and 
below (dashed lines) the energy at the Dirac point in a ID 
graphene superlattice. Red, green and blue lines correspond 
to the spatial period (L) being 10 nm, 15 nm and 25 nm, 
respectively, b, Similar quantities as in a. Red, green and 
blue lines correspond to the potential barrier width (w) being 
2.5 nm, 5 nm and 7.5 nm, respectively. 



Supplementary Discussion 6 : Fermi surfaces 

In ID and 2D graphene superlattices, the topology as 
well as the shape of the Fermi surface varies significantly 
with the Fermi energy (Supplementary Figs. 4 and 5). 
This variation gives rise to a dramatic variation in the 
species and the density of states of charge carriers as 
a function of the position of the Fermi level (Fig. 5 of 
paper) . 
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FIG. 8: (Supplementary Figure 4) Fermi surfaces of 
a ID graphene superlattice. a-f, Fermi surfaces of ID 
graphene superlattice with Uid = 0.3 eV, L = 10 nm and 
w — 5 nm plotted in the repeated zone scheme for different 
values of the Fermi energy (Ef) with respect to that at the 
Dirac point. Dashed lines are minizone boundaries. Red and 
blue lines are parts coming from the first and the second band 
above the Dirac point energy, respectively. 
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FIG. 9: (Supplementary Figure 5) Fermi surfaces of a 
2D rectangular graphene superlattice. a-f, Fermi sur- 
faces of a 2D rectangular graphene superlattice with U2D = 
0.3 eV, L x — L y — 10 nm and d = 5 nm plotted in the re- 
peated zone scheme for different values of the Fermi energy 
(.Ef) with respect to the Dirac point energy. Dashed lines 
are minizone boundaries. Red lines and blue line are parts 
coming from the first and the second band above the Dirac 
point energy, respectively. 



